function s=skew_matrix(a)
    s=zeros(3,3);
    s=[0    -a(3)   a(2);
       a(3)   0    -a(1);
       -a(2) a(1)     0];
end